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ABSTRACT 

We analyse the Tully-Fisher relation at moderate redshift from the point of view of the under- 
lying stellar populations, by comparing optical and NIR photometry with a phenomenological 
model that combines population synthesis with a simple prescription for chemical enrichment. 
The sample comprises 108 late-type galaxies extracted from the FORS Deep Field (FDF) and 
William Herschel Deep Field (WHDF) surveys at z <, 1 (median redshift z=0.45). A strong 
correlation is found between stellar mass and the parameters that describe the star formation 
history, with massive galaxies forming their populations early (zpoR ~ 3), with star forma- 
tion timescales, ri ^ 4Gyr; although with very efficient chemical enrichment timescales 
{t2 ^ 1 Gyr). In contrast, the stellar-to-dynamical mass ratio - which, in principle, would 
track the efficiency of feedback in the baryonic processes driving galaxy formation - does not 
correlate strongly with the model parameters. On the Tully-Fisher plane, no significant age 
segregation is found at fixed circular speed, whereas at fixed stellar-to-dynamical mass frac- 
tion, age splits the sample, with older galaxies having faster circular speeds at fixed Ms /A/dyn- 
Although our model does not introduce any prior constraint on dust reddening, we obtain a 
strong correlation between colour excess and stellar mass. 

Key words: galaxies: evolution - galaxies: formation - galaxies: stellar content - galaxies: 
fundamental parameters. 



1 INTRODUCTION 

Scaling relations among independent observational properties re- 
veal the presence of important mechanisms underlying the for- 
mation and evolution of galaxies, and the interplay between the 
baryons and the dark matter halos where galaxies reside. The strong 
correlation between luminosity an d circu lar speed in disc galax- 
ies (i.e. the Tully-Fisher relation, ITuIIv & Fisher 1977, hereafter 
TFR) has been used not onl y as an importan t rung in the cosmo- 
logical distance ladder (e.g. iGiovanelli et al.lll997,) . but also as a 
strong constraint on mod els of galaxy formation (e.g. iMathis et al.l 
uOOalDut ton et al. 2007). Its highest correlated version - the bary- 
onic TFR (McGaugh et al. 2000) - has eve n been used to te st the 
validity of standard Newtonian mechanics ( lMcGaughll2012n . The 
redshift evolution of the TFR allows u s to constrain the star forma- 
tion a nd assembly histories (see, e.g.. |Portinari & Sommer-LarsenI 
1200'^) . From the observational side, there are technical challenges 
that make the interpretation of the velocity field and the determi- 
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nation of M/L ratios rather complicated. After t he first results re - 
garding the redshift evolution of the TFR (e.g. IVogt et al.lll996r) . 
a nu mber of papers followed, claiming an evolution of the TFR 
slope jZiegler et al.ll2002l : iBohm et al.ll2004h . with low-mass disks 
being more l umino us in the the past. However, the analysis of 
IWeiner et al.l ( 120060 showed, instead, a significant brightnening of 
massive galaxies at high redshift, a result that could be recon- 
ciled if the star formation histories of massive galaxies have shorter 
timescales , as expected from m ulti-colour studies of disk galaxies 
(see, e.g.. iFerreras et al.ll2004l) . Constraining in a robust way the 
TFR is clearly a difficult task; the scatter of the relation is rather 
large for the amount of evolution observed, the outcome depends 
sensitively on the passband used, and sa mples at z^ 1 are inherently 
biase d towards the brigh test galaxie s ( Fernandez-Lorenzo et al.l 
|2009). Furthermore, iBohm & Ziegletj (l2007 b demonstrated that an 
evolution of the TFR scatter can mimic an evolution in slope. In 
addition, the complexity of disk kinematics is difficult to disentan- 
gle with single slit measurements, whereas 2D velocity fields from 
Integral Field Unit data suggest no evolution of slope, intercept 
or scatter of the TFR in disks with uniform kinematics, although 
the samples considered are rather sm all and also prone to biases 
dFlores et al.l2006l;|Puech et al.l2008h . 

Although semi-analytical models of galaxy formation predict 
a redshift evolution of the TFR in rough agreement with the ob- 
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servations, important pieces in the model ling of star formation are 
still missing (see, e.g. iTonini et al.ll201 11) . In this paper, we apply 
a simple phenomenological approach to gain insight on the TFR 
from the point of view of the underlying stellar populations. The 
more physically motivated version of the TFR, namely the cor- 
relation between circular speed and stellar mass, is an important 
indicator of disk galaxy formation. Comparisons with evolution 
models allow us to use stellar mass instead of luminosity, reduc- 
ing the highly variable changes in optical luminosity caused by 
the presence of young stars. The focus of this paper is to decon- 
struct the stellar mass T FR by exploring the age distribution of 
the stellar populations. In lFerreras et al.l (120041) we already showed 
that the stellar ages of disk galaxies at moderate redshift were a 
strong function of galaxy mass, possibly explained by a significant 
change in the star formation efficiency, with a strong increase at 
a circular speed Vmax > 140 km s^^. The mass-dependence of 
the stellar population properties of galaxies is a well-known phe- 
nomenon in the local universe (e.g. iKauffmann et a l. 2003;. Age 
indicators, like the strength of the 4000 A break ([Bruzual 1983), 
increase towards higher stellar masses, implying that high-mass 
galaxies formed sta rs with higher efficiency tha n low-mass ones 
in the cosmic past jPerez-Gonzalez et al.l 120081) . The connection 
between the mass of a galaxy and its stellar populations is also 
deeply intertwined with the evolution of the processes that drive 
galaxy formation over cosmic timescales. Overall, the main site 
of star formation shifts from high-mass galaxies at high redshift 
to successively lower-mass galaxies towards lower redshifts. This 
evolutionary trend is often referred to as downsizing (se e, e.g. 
ICowie et al.| [l996; Kodama et al. 20041: Icimatti et alj|2006l ). Sev- 
eral feedback processes have been proposed to explain such re- 
lationships, e.g. supernova feedback that regulate s star formation 
in low-mass galaxies (e.g. lOovemato et al.ll2009l) . or suppression 
of star formation b y active galactic nuclei in the center of high- 
mass galaxies (e.g. Khalatvan et alj|2008 ). It is still an ongoing de- 
bate whether the governing factor in such correlations is the stellar 
mass of a galaxy or rather the mass of its host dark matter halo. 
In the latter case, environment rather than mass would actually be 
the key factor that shapes the stellar populations. However, stud- 
ies base d on very large samples from the Sloan Digit al Sky Sur- 
vey (e.g. lvan den Bosch et alj2008l ; lRogers et al.l201(]f) present ev- 
idence that, at least at the present cosmic epoch, galaxy properties 
such as optical colour or morphological concentration are tightly 
correlated with stellar mass but only weakly with halo mass. 

In this paper, we revisit the evolution of disk galaxies by 
studying a sample at moderate redshift (z~0.5), with the aim of 
probing the evolution of the TFR and other scaling relations with 
respect to stellar ages. The outline of the paper is as follows: We 
present the photometric data from the FDF and WHDF samples 
(Section 2), and the phenomenological model to describe the star 
formation histories, including a comparison with more standard 
models of galactic chemical enrichment (Section 3), followed by 
a discussion of the main results (Section 4), along with our conclu- 
sions (Section 5). Throughout this paper, a standard ACDM cos- 
mology is assumed, with Q,m = 0.3 and Hq — 70kins~^Mpc~^. 



2 THE SAMPLE 

A sample of field disc galaxies is selecte d from the multi- band 
imaging of the FORS Deep Field (FD F, iHeidtetalJ 



William Herschel Deep Field (WHDF, iMetcalfe et al] 



20031) and 



200 ih 



veys. This sample was originally selected for the analysis of disc 



kinematics in the redshift range 0.1 < z < 1.0, using follow-up 
spectroscopic data t aken with the FORS camera at the Very La rge 
Telescope (see, e.g.. lB6hm e"tai]|2004 ISohm & Ziegleill2007h . In 
multi-object spectroscopy mode, a slit was placed along the appar- 
ent major axis of each target to extract a rotation curve, i.e. the 
rotation velocity as a function of radius. Hubble Space Telescope 
imaging with the Advanced Camera for Surveys (ACS) was used 
to determine parameters such as position angle, disc inclination 
or sca le length, with the use of the GALFIT package dPeng et al.l 
120021) . Taking into account all geometrical effects, such as inclina- 
tion, and all observational effects, such as seeing and slit width, the 
intrinsic maximum rotation velocity 1/max was determined for all 
galaxies with regular rotation curves. Neglecting objects with dis- 
turbed kinematics, solid-body rotation or a too low signal-to-noise 
ratio to probe into the regime of constant rotation velocity at large 
radii, Vmax could be determined for 124 galaxies, with a median 
uncertainty of A log Vmax = 0.087 dex. 

Photometry of our sample is available from the FDF and 
WHDF survey data. The FDF photometry has been acquired with 
the Very Large Telescope in the U-, B-, g-, R- and /-bands, and the 
New Technology Telescope in the J- and A's-bands. The WHDF 
photometry was carried out at the William Herschel Telescope 
in the U-, B-, R-, and /-bands and the Calar Alto 3.5m tele- 
scope in the H- and /i'-bands. Fixed aperture photometry was per- 
formed after convolving all frames to a common PSF with FWHM 
1 (1.5)arcsec for FDF (WHDF) galaxies. The aperture diameter 
used for the photometry is 2 (3) arcsec for FDF (WHDF) sources. 
The fi nal magnitudes were corrected for Galactic extinction fol- 
lowing Ic^eUTetal] ( Il989l) . adopting a reddening E{B — V) = 
0.018 mag and E{B — V) = 0.030 mag to wards the positions o f 
the FDF and WHDF, respectively (based on ISchlegel et ai]|l998h . 
A number of galaxies with available Vmax values were not ob- 
served in all filters and thus rejected, resulting in a final sample of 
108 disks (73 in FDF and 35 in WHDF), between z = 0.07 and 
z = 0.97 at a median redshift of (z) = 0.45. 

Star formation rates (SFRs) were estimated from [O II] equiv- 
alent widths, available for 65 galaxies, following iKennicutll 119921) 
(other emission lines, such as H/3 or [O III] are sometimes used for 
the determination of the rotation curves, depending on the wave- 
length range probed by a given spectrum). Note that Ha is not avail- 
able within the observed wavelength range to derive SFRs given the 
redshift distribution of the sample. 

Figure [T] shows a comparison of some of the available 
colours with respect to two sets of synthetic models from 
ISruzual & CharlotI ( l2003h : simple stellar populations (i.e. single 
age and metallicity, left) or composite models with an exponentially 
decaying star formation history, at fixed metallicity (i.e. r-models, 
right). Different lines correspond to choices of stellar age in SSP 
models (given by the formation redshift), or exponential timescale 
and metallicity for the r models, as labelled. The combination of 
colours across a wide range of wavelengths allows us to constrain 
the ages and metallicities: the optical colours (top) are most sensi- 
tive to age, whereas the NIR colours have an increased dependence 
with respect to metallicity. All these models are dustless, except for 
the thin solid line in the r models, which co rresponds to an intri nsic 
reddening E{B - V) = 0.2 (following the lCardelliet al.ll989l. ex- 
tinction law). The comparison shows that the reddening values stay 
below E{B — V) ;$ 0.5 mag. The figure also illustrates the need 
to choose composite models to explain the multi-band photometric 
data, as expected from the complex star formation histories of disc 
galaxies. 
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Figure 1. Comparison of the photometry with simple stellar populations (left) and models with exponentially decaying star formation histories at fixed 
metallicity (right), all based on the Brazual & Chariot (2003) population synthesis models. The thin solid line in the EXP panels corresponds to a colour 
excess of E(B — V) = 0.2 mag for HTf = 5 Gyr model, at solar metallicity. Galaxies in the PDF (WHDF) are shown as sohd (open) dots, where a small 
offset (always below 0.1 mag) is applied to the WHDF data to correct for the different response of the filters used (only for display purposes in this figure). 



3 MODELLING STAR FORMATION HISTORIES 

We apply a phenomenological one-zone model that describes the 
star formation history (SFH) and chemical enrichment of a galaxy 
with three parameters. The model is constrained by the broadband 
photometric data presented above. A simple model allows us to 
probe a large volume of parameter space, deriving more consistent 
constraints on the mechanisms contributing to the build-up of the 
stellar populations. The star formation history is assumed to start at 
an epoch given by a formation redshift (zfor), with a star forma- 
tion rate modelled either by a delayed exponential: 



or by a standard exponential (i.e. a r-model): 
© 2013 RAS, MNRAS 000,[T]] 



(1) 



(2) 



and a metallicity trend driven by the growth in stellar mass, al- 
though with an independent timescale, to take into account the ef- 
fects of both gas outflows and star formation efficiency: 



Z{t) = Zi+Z2\l- exp 



-[t-t(zFOR)l/T2 



(3) 



The extrema in metallicity are Zi = Zq/10 and Z2 — 2Zq. In 
this model, n gives an indication of the timescale for the formation 
of the stellar component, whereas T2 drives the metal enrichment. 
These two timescales are considered independent parameters in the 
analysis, hence we allow for scenarios with a decoupling between 
star formation, gas infall and outflows. As an example, Fig.|2]shows 
the relationship between star formation rate (assuming a delayed 
exponential law) and metallicity for three models. Note that the dot- 
ted line, corresponding to a balance between enrichment and infall, 
gives a metallicity distribution similar to the closed box model, pro- 
ducing an excess of stars at low metallicity. This so-called G-dwarf 
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Figure 2. The model adopted to track the star formation histories of disc 
galaxies includes two different timescales: ri tracks the star formation rate 
- which evolves as a delayed exponential. T2 controls the emichment rate. 
This figure shows the difference in the time evolution of the star formation 
rate and the metallicity ([m/H]) for three cases, as labelled. An additional 
parameter (formation redshift) controls the average age of the distribution. 



problem (see, e.g..' Pagelll 19971) - is solved in this framework by a 
fast enrichment (short r2) along with an extended period of infall 
(long ri). In Fig.[3l we compare the age and metallicity distribu- 
tions given by this model with a more detailed a nalysis based on 
standard prescriptions of c hemical enrichment (see lFerreras & SilM 
l2001uFerreras et al.l l20Q4, for details). Contours of average stellar 
age (black lines) and metallicity (grey lines) are shown with respect 
to star formation efficiency (i/) and outflow fraction (/3) for the en- 
richment model (Jeft) or with respect to the two timescales (ri and 
T2) for the phenomenological model used here (right). The forma- 
tion epoch is fixed in both cases to zpor = 5. A mapping can be 
established between star formation efficiency (v) and our ri ; and 
also between outflow fraction (/3) and our r^. We emphasize that 
the model presented here is not only simpler, but it also gives a 
clearer description of the formation of the stellar populations. 



4 RESULTS 

The models presented in the previous section are explored over a 
large range of parameters. We also include the effects of dust as 
a homogeneous screen. We describe this component as an addi- 
tional fourthfree_£arametet^hoosing the colour excess, E{B~V), 
from the lCardelli et alj ( 119890 extinction law. In order to obtain the 
most robust estimates of the best fit and uncertainty, we opt for 
a comprehensive search of parameter space. This is doable given 
the small number of parameters describing the star formation his- 
tory. Furthermore, a parameter search algorithm has a very simple 
parallelisation scheme, so that we can run it efficiently on many 
processors. For a choice of the four parameters, the resulting star 
formation history i s com bined with the synthetic populations of 
ISruzual & CharlotI (l2003h to derive the six (five) broadband pho- 
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Figure 3. Comparison of the model used in this paper with a more generic 
chemical enrichment model that includes parameters to describe the star 
for mation efficiency (1/) an d the fraction of gas ejected in outflows (J3) (see 
e.g. iFeneras & SilklbOOlL for details). In order to simplify the compari- 
son, we assume the same formation epoch in both cases (zpoR = 5). The 
panel on the left shows contours of mass-weighted average age (black) and 
metallicity (grey), with respect to star formation efficiency and gas outflow 
fraction. The black (grey) labels indicate the extrema in age (metallicity). 
The panel on the right shows similar contours for the model used here, ex- 
ploring a range of values of star formation timescale (ri) and enrichment 
timescale (t2). 



tometric colours available from the FDF (WHDF) data, generating 
a likelihood (via a standard x^ statistic), that defines a probability 
distribution function for the three parameters under consideration. 
Out of the two families of models - corresponding to either a de- 
layed exponential, or a r model - we select for each galaxy the 
one that gives the lowest value of the x^- Although discriminat- 
ing between these two prescriptions for the star formation history 
is beyond the capabilities of broadband photometry, we note that 
there is a slight preference towards the delayed exponential model, 
especially at low stellar mass. We derive from this likelihood the 
average mass-weighted age and metallicity, along with the uncer- 
tainties, quoted throughout this paper at the 69% level. For each 
galaxy, we run two grids of 32^ x 16 models (one for each choice 
of star formation rate), uniformly scaled over the following range 
of the parameters: 

1 <zfor < 5 
-ls;log(ri/Gyr)s;+l 
-l^log(T2/Gyr)s;+l 

0<£(B- V) s;0.5 

In addition, we impose a prior on the average stellar metal- 
licity, using the local mass-metallicity relation of lOallazzi et al.l 
(2005). We note that this prior is rather mild, given the large 
scatter of the distribution at fixed mass (A[Z/H] ~ Idex at 
log Ms/Mq = 10). The reason for imposing this prior is based on 
the lack of a large number of independent metallicities in the syn- 
thetic models, along with the fact that broadband photometry alone 
- only extending out to rest-frame data in the H band, given the red- 
shift range - cannot impose strong constraints on the metallicities 
of unresolved stellar populations. This prior helps to ensure a cor- 
rect age-metallicity discrimination in the analysis. The error bars 
on the estimates of stellar metallicity are around A[Z/H] ~ 0.3 dex 
and the average age uncertainties A(Age)/(Age) = 0.4. A com- 
parison between the ages derived with and without the prior on 
metallicity gives a difference (Age) prior - (Age) no prior = 
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Figure 4. Recovery of average (top) and RMS of the age distribution (bot- 
tom) from 200 simulations spanning a similar range of redshift and SNR 
as the original data. The vertical axis in each panel represents input mi- 
nus reco vered value. The blac k (grey) data points use the metallicity prior 
from the lGallazzi et al.l 420051) relation, assuming a stellar mass of 10^^ M© 
(10^ Mq). Error bars are the 69 % confidence levels from the analysis, and 
the histogram of the distribution is shown in the rightmost panels. The val- 
ues of the difference in average age and RMS of each case is indicated with 
the same colour coding. 



Figure 5. Recovery of model parameters from 200 simulations spanning a 
similar range of redshift and SNR as the original data. In each case A is de- 
fined as the difference between the input and the recovered parameter. The 
black (grey) histograms con'espond t o the application of th e prior regarding 
the mass-metallicity relationship of lGallazzi et al.l 4200511 assuming galax- 
ies with stellar mass of 10^ Mq (10^^ Mq), to illustrate the effect of such 
a prior on the recovery of the model parameters (see text for details). The 
numbers in each panel show the RMS of the distributions for each choice 
of stellar mass, with the same colour coding. 



0.14 ± 0.61 Gyr, without any systematic effect with respect to 
galaxy mass. 

In order to confirm this point, we performed 200 simulations 
with the same redshift and signal-to-noise distribution as the origi- 
nal data, comparing the age and RMS of the age distribution. Fig.|4] 
shows the difference between the input and the recovered values for 
the average age (top) and the RMS of the age distribution (bottom). 
The simulations assume two different cases of stellar mass for the 
application of the prior: 10^" Mq (black) and 10^ M© (grey). In 
either case, the recovered values do not show any systematic off- 
set, with an accuracy of ~ 0.9 Gyr in average age and --^0.3 Gyr in 
the RMS of the age distribution. In Fig. |5]a similar comparison is 
made, in this case between the model parameters. The numbers in 
each panel give the RMS of the distributions of the (input — out- 
put) values of the parameters for the simulations. Note that there is 
no systematic change in the values of the parameters, except for T2. 
In this case, as expected, the prior biases the values of metallicity 
towards lower metallicity for the lO" Mq (grey) case, and towards 
higher metallicity at lO" Mq (black)B 



4.1 Model parameters 

Hereafter, we present the results for the FDFAVHDF data. Fig.|6] 
shows the distribution of the two formation timescales. The sample 



We note that the prior for the simulations is blindly applied in all galaxies, 
irrespective of their input value of metallicity, hence representing a worst 
case scenario. 
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Figure 6. The model timescales t\ (star formation) and T2 (enrichment) 
are shown for our sample of FDF+WHDF disc galaxies, split with respect 
to stellar mass, as labelled. A typical 1 a en'or bar is shown. 



is split with respect to stellar mass, at the position of the median 
(~ 10^'^Mq). For reference, the median uncertainty (at the 68% 
level) is shown in each panel as a horizontal error bar. Note the 
difference between the distribution of low- and high-mass discs, 
especially in the enrichment timescale, T2. The majority of low- 
mass discs (grey histograms) are fitted by longer timescales both 
in star formation (t\) and enrichment (r2). A significant fraction of 
the massive galaxies are better explained by very short enrichment 
timescales (t2 ~ 1 Gyr). Low values of both n and T2 would repre- 
sent the efficient and fast buildu p of the stellar populati ons expected 
in early-type galaxies (see, e.g.. lde la Rosa et al.l201ll) . In contrast, 
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Figure 7. The constraints on the parameters of the models are shown as a function of total mass (Jeff), stellar mass (centre) and stellar-to-total mass fraction 
(right). Individual values are shown as dots, whereas the black crosses represent the median values after binning (with equal number of galaxies per bin), 
including the RMS scatter per bin as an en'or bai'. The red solid (blue open) dots coiTespond to ^ ^ 0.45 (< 0.45) (i.e. split at the median redshift of the 
sample). 



massive discs have more extended star formation timescales (re- 
flecting a lower star formation efficiency or an extended infall), and 
short enrichment timescales (possibly caused by a lack of metal- 
rich outflows). The difference in the error bars for n and T2 is 
caused by the higher accuracy on the age estimates, with respect 
to metallicity. 

Fig. |7] shows the best fit parameters as a function of physical 
properties: dynamical mass (left), stellar mass (centre), or stellar- 
to-dynamical mass fraction (right). To avoid model-dependent ex- 
trapolations for the dark matter distribution, both dynamical and 
stellar masses are quoted within a fixed aperture of 3 exponen- 
tial scale lengths (roughly where all galaxies in the sample reach 
a flat circular velocity). The scale le ngths are derived from surface 
brightness fits to the ACS images JBohm & Ziegleiil2007l) . For a 
clearer interpretation of the trends, the individual data points are 
arranged into four bins. The binning is done at constant number of 
galaxies per bin, and the median and dispersion (RMS) are shown 
as crosses, and error bars, respectively. The individual data points 
are also coded with respect to redshift (red solid: z Js 0.45; open 



blue: z < 0.45). Note that the parameters are correlated with stel- 
lar mass, with a clear downsizing trend with zpor. The correlations 
with respect to Ms/Mdyn could be explained with the assumption 
that this fraction is an indicator of efficiency or gas outflows, hence, 
at large stellar masses, one should expect a fast, efficient buildup 
of metallicity (i.e. a short T2]. Although such trend is present in 
the data, its correlation is much weaker than with respect to stel- 
lar mass. Regarding redshift, the strongest trend is with star forma- 
tion/infall timescale (ri, middle panels), reflecting a more extended 
distribution of ages in the sample at low redshift, mostly due to the 
wider range of ages because of lookback time. 

The trends of average age and metallicity with stellar mass are 
shown in Fig.[8l split into three redshift bins, that correspond to uni- 
form steps of ~ 2Gyr in lookback times between bins. Error bars 
are shown at the 68% confidence level. A linear fit between age and 
mass is applied to the high redshift bin (top left panel), and shifted 
by the lookback time in the other two bins, for reference. Within 
error bars, the slopes are consistent, with the only change being 
an offset amounting to the change in lookback time. Although this 
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Figure 8. Distribution of ages (Jeft) and metallicities (right) of the 
FDF+WHDF disc galaxies, split with respect to redshift in three panels. 
The straight Hne on the top-left panel fits the data points in the high red- 
shift subsample. This line is copied over, as a dashed hne, in the other two 
redshift bins (middle- and bottom-left panels). We also show the effect of 
shifting this fit by the lookback time difference, as solid lines in these two 
panels. On the right, the solid and dashed line s mark the average a nd 1 a 
unceitainties of the mass-metalhcity relation of lGallazzi et alj 120051) . 



could be interpreted as a lack of "dynamical downsizing", the sam- 
ple is not large enough to constrain a change in slope of the cor- 
relation between age and stellar mass. As regards to metallicity, 
there is a trend towards the metal rich envelope of the local mass- 
metallicity relation. Note the present analysis, based on broadband 
photometry alone, mainly constrains stellar ages, whereas estimates 
of metallicity carry a large error bar (typically A[Z/H] ~ 0.3 dex). 
Even though the local mass-metallicity relation is used as a mild 
prior in the analysis, we stress that the minimum values of x^ are 
not significantly affected by the prior, with ~ 90% of the sample 



having a reduced value x'. 



2 < 



Figure 9. Relationship between maximum rotation velocity and stellar mass 
(left) or stellar-to-dynamical mass fraction (right). The sample is split at the 
median with respect to stellar average age. Top and bottom panels show 
subsamples at high and low redshift, respectively (also split at the median 
of the distribution). The sohd hne in each panel is the best fit to the data 
within each redshift bin. Typical error bars are shown. The dashed line in the 
left panels show t he local rel ationship from the 7-band determined stellar 
masses I Bell & d e JonalOOlh . The dotted line is the z~l relationship from 
jMilleretal.l20lJl 



a faster Vmax for disc galaxies with a lower stellar mass fraction, 
as expected, since the stellar-to -total mass ratio d ecreases with in- 
creasing galaxy mass (see, e.g.. lMoster et al.ll201 0). At fixed mass 
fraction, the younger galaxies have slower rotation speeds. One 
can assume that the stellar-to-dynamical mass fraction controls the 
level of "baryon feedback" - mainly related to star formation in this 
sample. A tentative interpretation of the segregation with respect to 
age, at fixed AIs/M^yn, would be that older galaxies are preferen- 
tially assembled in higher density haloes, therefore having a higher 
circular speed. 



4.2 The evolution of the Tiilly-Fisher relation 

Fig. [9] shows the properties of the under lying stellar populatio ns on 
a variation of the TuUy-Fisher diagram (ITullv & Fisheijl 19771) . The 
sample is split at the median redshift (top/bottom panels), and in 
each panel, black solid (grey open) circles represent disc galaxies 
with an average age older (younger) than the mean value within 
each redshift bin, as labelled. Typical (la) error bars are shown 
as reference. The panels on the left show a Vmax vs stellar mass 
TuUy-Fisher relation. The usual downsizing relation between stel- 
lar mass and age is found. However, at fixed stellar mass, no signif- 
icant segregation with respect to Vmax is evident. Old and young 
discs show a similar amount of scatter with respect to the best lin- 
ear fit (solid line). The slope and intercept a re indistinguishable (as 
also found by, e.g.'Miller et al.ll201ll [20121 , out to Z--1.7) with re- 
spect to the redshift subsample chosen (Ms oc Vj^ax)- The pan- 
els on the right show the maximum circular speed with respect 
to stellar-to-dynamical mass fraction, with a clear trend towards 



4.3 The contribution of dust 

Fig. [To] shows the trend of dust reddening with respect to stellar 
mass, average metallicity or inclination. The sample is split at the 
median value of the stellar mass (^ 10^" Mq), with the red solid 
(blue open) dots corresponding to stellar masses above (below) 
this cut. The big solid dots with error bars in the left and middle 
panels give the median and RMS scatter of the complete sample, 
binned at fi xed number of galaxies per bin. A prescription from 
ITuIIv & Fo uque ( 1985) is included for reference, with / — 0.2 and 
T — 0.5 (solid line). A strong trend is found with respect to mass 
and metallicity, in good agreement with the observations of local 
disks jTullv et al.lll998l) . More subtle is the excess of more attenu- 
ated galaxies at higher inclinations for the more massive galaxies. 
Notice that our modelling does not impose any prior constraint on 
the colour excess, con sidering all values of E (B — V) in the anal- 
ysis. The formalism of ITuUv & Fougug ( Il985n is applicable for the 
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Figure 10. Model constraints on the colour excess, with respect to stellar mass, metallicity or inclination. The red solid (blue open) dots correspond to stellar 
masses log Ms/Mq > 10 (< 10). The big solid dots with eiTor bars in the left and middle panels give the median and RMS scatter of the complete sample, 
binned at fixed number of galaxies per bin. The rightmost panel includes a prescription from lTullv & Fouqug ( Il985b . with / = 0.2 and r = 0.5 (solid line). 



massive disc subsample, whereas lower mass discs are not affected 
so strongly with respect to disc inclination. 



4.4 Star formation rates (SFRs) 

Ongoing star formation rates can be derived from emission line 
luminosities. In this paper, we use SFRs from [O II], using the con- 
version from lKennicuttI l ll992h . Fig.[TT]presents the SFRs (top) and 
specific SFRs (SSFR, defined as the SFR per unit stellar mass, bot- 
tom) with respect to either stellar mass {left) or average age (right). 
The sample is spUt with respect to redshift at the median value of 
the subsample comprising 65 galaxies where SFRs could be de- 
termined (zm ~ 0.60), with solid (open) circles corresponding 
to low (high) redshift. Only a subset of the full sample is pre- 
sented here because the derivation of SFRs require estimates of 
[On] luminosities. Characteristic error bars, at the la level are 
shown for reference. A strong correlation is found between S SFR 
and age. In the bottom-left panel, the rates measured in our sample 
are compared with the gene ral trend of star- forming galaxies from 
iBrinchmann & EllisI feOOd) at 2 ~ 0.6 (solid line) and z ~ 0.9 
(dashed line). 



(see, e.g. iFerreras & SilkiBOOlh . however, our massive discs have 
longer SF timescales, resulting in an extended distribution of stel- 
lar ages. On the TuUy-Fisher plane, age does not introduce a seg- 
regation at fixed circular speed, suggesting mild evolution of the 
Tully-Fisher relation with redshift. 

We also present the model fits with respect to dynamical mass 
(measured within three scale lengths) and find a similar trend as 
with respect to stellar mass. However, the ratio of the two - an 
estimator of dark matter fraction, or feedback efficiency - seem 
not to correlate so strongly with the model parameters, reinforc- 
ing the idea that mass is the strongest, first order driver of the star 
formation histories of galaxies. Nevertheless, on the Tully-Fisher 
diagram (Fig.|9]l, at fixed Ms/Mdyn, older galaxies have faster cir- 
cular speeds, perhaps reflecting the fact that older discs would have 
formed in higher density halos, hence the higher Vmax. The dif- 
ference is rather small, and a more detailed analysis - based on 
spectroscopic data with high SNR and accurate flux calibration - is 
needed to confirm this trend. 

Our model includes dust as a free parameter, treated as a sim- 
ple dust screen. Our findings agree independently with the formal- 
ism of iTullv & Fouqug ( Il985h for the most massive discs (stellar 
mass ^ 10^^ Mq), with lower mass discs featuring no significant 
trend between colour excess and inclination. 



5 CONCLUSIONS 

By comparing optical to NIR photometry of a sample of disc galax- 
ies at intermediate redshift, with a simple set of phenomenological 
one-zone models including chemical enrichment, we are able to ex- 
plore the formation of these systems and the evolution of the Tully- 
Fisher relation. We find a strong correlation between the parameters 
that describe the star formation history and stellar mass. The forma- 
tion epoch and star formation (SF) timescale (ri) are especially cor- 
related with stellar mass, showing the usual downsizing trend, with 
the most massive discs forming around z ~ 3, with SF timescale 
ri ~ 4Gyr. The enrichment timescale (T2) decreases with increas- 
ing mass, reflecting a very efficient process for the build-up of 
metallicity in the most massive galaxies (t2 ;$ 1 Gyr). These are the 
typical timescales for the enrichment of massive early-type galaxies 
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Figure 11. Correlation of star formation rate (SFR, top) and specific star 
formation rate (SSFR, bottom) witli respect to stellar mass (left) or average 
age of the stellar populations (right). The data are split at the median red- 
shift of the subsample where SFR information is available (red solid dots: 
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